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Superelasticity  (SE)  in  bulk  materials  is  known  to  originate  from  the  structure-changing  martensitic  transi¬ 
tion  which  provides  a  volumetric  thermodynamic  driving  force  for  shape  recovery.  On  the  other  hand, 
structure-invariant  deformation  processes,  such  as  twinning  and  dislocation  slip,  which  result  in  plastic  defor¬ 
mation,  cannot  provide  the  driving  force  for  shape  recovery.  We  use  molecular-dynamics  simulations  to  show 
that  some  bcc  metal  nanowires  exhibit  SE  by  a  “reversible”  twinning  mechanism,  in  contrast  to  the  above 
conventional  point  of  view.  We  show  that  this  reversible  twinning  is  driven  by  the  surface  energy  change 
between  the  twinned  and  detwinned  state.  In  view  of  similar  recent  findings  in  fee  nanowires,  we  suggest  that 
SE  is  a  general  phenomenon  in  cubic  nanowires  and  that  the  driving  force  for  the  shape  recovery  arises  from 
minimizing  the  surface  energy.  Furthermore,  we  find  that  SE  in  bcc  nanowires  is  unique  in  several  respects: 
hrst,  the  (111)/{11 2}  stacking  fault  generated  by  partial  dislocation  is  always  preferred  over  (111}/{110}  and 
(1 1 1)/{123}  full  dislocation  slip.  The  occurrence  of  (111)/{112}  twin  or  full  dislocation  slip  in  bcc  nanowires 
depends  on  the  competition  between  the  emission  of  subsequent  partial  dislocations  in  adjacent  {112}  planes 
and  the  emission  of  partial  dislocations  in  the  same  plane.  Second,  compared  to  their  fee  counterparts,  bcc 
nanowires  have  a  higher  energy  barrier  for  the  nucleation  of  twins,  but  a  lower  energy  barrier  for  twin 
migration.  This  results  in  certain  unique  characteristics  of  SE  in  bcc  nanowires,  such  as  low  energy  dissipation 
and  low  strain  hardening.  Third,  certain  refractory  bcc  nanowires,  such  as  W  and  Mo,  can  show  SE  at  very  high 
temperatures,  which  are  higher  than  almost  all  of  the  reported  high-temperature  shape  memory  alloys.  Our 
work  provides  a  deeper  understanding  of  superelasticity  in  nanowires  and  refractory  bcc  nanowires  are  poten¬ 
tial  candidates  for  applications  in  nanoelectromechanical  systems  operating  over  a  wide  temperature  range. 

DOI:  10.1103/PhysRevB. 82.205435  PACS  number(s):  62.25. -g,  62.23.Hj,  61.72.Mm,  68.35.Md 


I.  INTRODUCTION 

Superelasticity  (SE)  is  a  characteristic  property  of  shape- 
memory  alloys  (SMAs)  which  have  been  widely  studied  over 
several  decades  due  to  their  applications  as  sensors,  actua¬ 
tors,  and  microeletromechanical  devices.1^3  In  general,  su¬ 
perelasticity  is  based  on  the  phase  stability  of  thermoelastic 
martensite:  a  stress-induced  martensitic  transformation  dur¬ 
ing  loading  and  the  reverse-martensitic  transformation  dur¬ 
ing  unloading.4  In  contrast,  deformation  of  bulk  crystals  is 
often  controlled  by  the  motion  of  dislocations  or  deformation 
twins,  processes  that  result  in  plastic  deformation  and  thus 
cannot  usually  recover  after  unloading. 

However,  nanosized  materials  possess  the  feature  that  the 
surface  energy  (negligible  for  bulk  materials)  becomes  im¬ 
portant  and  this  can  lead  to  novel  properties  not  expected  in 
bulk  materials. 5-9  Recently,  molecular-dynamics  (MD)  simu¬ 
lations  have  shown  that  a  class  of  fee  metals 10-1 3  in  the  form 
of  nanowires  can  exhibit  SE  via  a  twinning  mechanism.  For 
example,  in  single-crystalline  Cu,  Ag,  and  Ni  metallic  nano¬ 
wires,  the  initially  created  ( 100) / {100}  nanowires  can  spon¬ 
taneously  reorient  to  a  defect-free  (110) /{111}  structure 
above  a  critical  temperature.10-12-15  Interestingly,  under 
uniaxial  tension,  the  reoriented  ( 1 10)/{1 11}  wire  can  exhibit 
SE  by  transforming  to  the  initial  (100)/{100}  configuration 
via  movement  of  twin  boundaries.  u~15  This  mechanism  for 
SE  has  not  been  previously  considered  and  is  quite  different 


from  that  associated  with  bulk  SMAs.1  It  has  been  well  ac¬ 
cepted  that  the  large  surface  energy  difference  between  the 
twinned  and  detwinned  configurations  as  well  as  the  move¬ 
ment  of  twins  are  the  two  major  factors  that  dominate  this 
unique  behavior.  u)~12-14  The  minimization  of  surface  energy 
for  a  nanowire  provides  the  thermodynamic  driving  force 
whereas  the  reversible  twin  boundary  motion  provides  a 
crystallographic  pathway  for  the  complementary  reorienta¬ 
tion.  However,  this  phenomenon  has  mainly  focused  on  fee 
nanowires  to  date  and  it  is  unclear  whether  this  effect  is 
common  to  nanowires  possessing  other  crystal  structures,  in 
particular  nonclose  packed  lattices. 

The  purpose  of  this  work  is  to  use  the  MD  simulation 
approach  to  explore  superelasticity  in  bcc  nanowires  via  a 
similar  twinning  mechanism  as  in  fee  nanowires.  Our  results 
show  that  some  bcc  nanowires  also  deform  superelastically 
via  a  mechanism  that  involves  minimization  of  surface  en¬ 
ergy  and  twin  boundary  motion.  Our  work  emphasizes  that 
superelasticity  is  a  general  phenomenon  in  cubic  metal  nano¬ 
wires  where  both  the  large  surface  energy  difference  and 
reversible  twinning  deformation  play  key  roles.  Moreover, 
we  find  that  in  contrast  to  fee  nanowires,  the  SE  of  bcc 
nanowires  is  unique  in  several  respects:  (i)  the  (111)/{112} 
stacking  fault  generated  by  the  |(111)  partial  dislocation  is 
always  preferred  to  the  (111)/{110}  and  (11 1) /{ 123}  full  dis¬ 
location  slip.  The  occurrence  of  the  (111)/{1 12}  twin  or  full 
dislocation  slip  in  bcc  nanowires  depends  on  the  competition 
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between  the  emission  of  subsequent  §(111)  partial  disloca¬ 
tions  in  adjacent  {112}  planes  and  the  emission  of  §(111) 
partial  dislocations  in  the  same  {112}  plane,  (ii)  bcc  nano¬ 
wires  face  a  higher  energy  barrier  for  nucleating  twins  but  a 
lower  barrier  for  twin  migration,  compared  to  their  fee  coun¬ 
terparts.  This  results  in  unique  characteristics  related  to  SE  in 
bcc  nanowires,  such  as  low  energy  dissipation  and  low  strain 
hardening,  (iii)  Certain  refractory  bcc  nanowires,  such  as  W 
and  Mo,  show  SE  at  very  high  temperatures,  which  are 
higher  than  almost  all  of  the  reported  high-temperature 
SMAs. 16,17  Hence,  our  work  has  the  potential  to  impact  ap¬ 
plications  of  materials  showing  superelasticity  in  nanoelec¬ 
tromechanical  systems  operating  over  a  wide  temperature 
range. 

The  outline  of  the  paper  is  as  follows.  In  Sec.  II,  we 
present  the  simulation  method,  whereas  Sec.  Ill  focuses  on 
the  simulation  results.  In  Sec.  IV,  we  analyze  the  driving 
force  for  superelasticity  in  bcc  nanowires,  and  consider  the 
competition  between  twinning  and  full-dislocation  slip  in 
these  nanowires.  Superelasticity  in  fee  and  bcc  nanowires  is 
then  compared  and  we  conclude  with  Sec.  V. 

II.  SIMULATION  METHOD 

In  our  MD  simulations  we  focus  on  nanowires  of  iron,18 
vanadium,19  tungsten,20  and  molybdenum20  using  embedded- 
atom  method  (EAM)  (Ref.  21)  potentials  which  reproduce 
relatively  good  thermomechanical  properties  of  the  above 
bcc  element  metals,  such  as  the  cohesive  energy,  elastic  con¬ 
stants,  and  the  low-index  surface  energies  (especially  the  sur¬ 
face  energy  difference  between  the  low-index  surfaces, 
which  play  a  crucial  role  in  the  present  calculations),  etc. 18-20 
However,  one  should  not  expect  that  such  semiempirical  po¬ 
tentials  can  reproduce  all  the  properties  of  the  metals.  The 
initial  bcc  nanowires  were  first  created  along  the  x 
-[100],y-[011],z-[0ll]  directions  at  0  K.  The  wires  have 
dimensions  of  50a  X  5 \2a  X  5  \  2a,  where  a  is  the  lattice  con¬ 
stant  at  0  K  (0.31472  nm,  0.31652  nm,  0.28553  nm,  and 
0.30299  nm  for  Mo,  W,  Fe,  and  V,  respectively).  We  then 
relaxed  the  initial  nanowires  at  300  K  for  almost  1000  pico¬ 
seconds  to  obtain  an  equilibrium  state  by  using  the  Nose- 
Hoover  thermostat.22-2'  The  relaxed  wires  were  then  loaded 
in  the  x  direction  with  atoms  at  the  ends  fixed  to  play  the  role 
of  loading  grips.  The  tensile  strain  was  applied  by  a  ramp 
velocity  profile  to  give  the  minimum  strain  at  one  fixed  end 
and  maximum  at  the  other.  Unloading  was  performed  in  a 
similar  manner  by  applying  a  reverse  ramp  velocity  that  was 
terminated  when  the  stress  was  reduced  to  zero.  The  strain 
rate  during  loading/unloading  was  107  s_1  and  the  MD  cal¬ 
culations  were  carried  out  using  the  LAMMPS  code24  with 
the  atomic  configurations  displayed  by  ATOMEYE.25 

III.  RESULTS 

A.  Superelasticity  in  Mo,  W,  and  Fe  nanowires 

The  mechanical  behavior  of  a  Mo  nanowire  was  first 
tested  under  unixial  tension.  Figure  1(a)  shows  the  stress- 
strain  curve  of  the  Mo  nanowire  upon  loading  and  unloading 
at  300  K.  We  see  that  in  the  loading  stage,  after  a  short 
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FIG.  1.  (Color  online)  Illustration  of  superelasticity  at  300  K 
in  a  Mo  nanowire  with  dimensions  of  15.74  nmX2.23  nm 
X  2.23  nm.  (a)  Stress-strain  curve  of  the  Mo  nanowire  during  the 
loading  and  unloading  process.  The  recoverable  strain  can  exceed 
40%.  (b)  The  corresponding  atomic  configurations  due  to  mechani¬ 
cal  deformation  (only  two  adjacent  lattice  planes  are  shown  for  each 
snapshot).  Upon  loading,  the  [100]-axis  configuration  transforms  to 
the  [011]-axis  configuration  with  the  aid  of  twinning.  The  lattice 
unit  cell  rotates  by  90°  to  produce  the  crystal  reorientation.  Subse¬ 
quently,  upon  unloading  the  [011]-axis  nano  wire  reverts  back 
through  reversible  twin  boundary  motion.  Atoms  are  colored  ac¬ 
cording  to  the  values  of  their  potential  energy  in  the  unit  of  eV. 

elastic  deformation  to  the  yield  point  of  16.2  GPa,  the  stress 
drops  dramatically  to  a  relatively  low  value.  Further  loading 
causes  an  extensive  plateau  stage  with  stress  remaining  at 
approximately  2.2  GPa.  Upon  unloading,  the  plastic  strain 
can  almost  fully  recover  along  the  same  pathway  as  the  load¬ 
ing  process  to  exhibit  classic  superelastic  behavior  with  a 
total  recoverable  strain  exceeding  40%. 

By  observing  the  corresponding  atomic  configurations 
during  loading/unloading,  we  find  that  twinning/detwinning 
is  the  underlying  mechanism  responsible  for  the  reversible 
strain.  As  shown  in  Fig.  1(b),  upon  tensile  loading  the  newly 
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formed  (111)/{112}  single  twin  divides  the  whole  nanowire 
into  three  domains.  As  strain  increases,  the  twin  boundaries 
further  propagate  toward  the  ends  and  the  whole  wire  trans¬ 
forms  from  the  initial  orientation  (x-[100],y-[011],z 
—[Oil])  to  the  new  configuration  (x-[011],y-[100],z 
-[Oil])  with  a  90°  rotation  of  the  unit  cell  in  the  xy  plane. 
For  convenience,  we  hereafter  refer  to  the  initial  and  reori¬ 
ented  configurations  of  the  wire  as  [100]  and  [Oil],  respec¬ 
tively.  When  the  system  is  unloaded,  detwinning  follows  the 
reverse  path  to  allow  the  wire  to  recover  its  original  shape. 
The  presently  observed  SE  due  to  twinning  is  similar  to  that 
shown  in  fee  nanowires. 10-14 

We  further  investigated  the  mechanical  behavior  of  W  and 
Fe  nanowires  under  uniaxial  tension.  Figures  2(a)  and  2(b) 
show  the  stress-strain  curves  of  W  and  Fe  wires  at  300  K, 
respectively.  From  Figs.  2(a)  and  2(b),  we  clearly  see  that  the 
SE  also  exists  in  W  and  Fe  nanowires.  Atomic  configurations 
(not  shown  here)  corroborate  that  the  SE  of  W  and  Fe  nano¬ 
wires  is  also  mediated  by  a  twinning/detwinning  mechanism, 
as  for  Mo  nanowires.  One  interesting  observation  is  that  the 
stress-strain  curves  for  Mo,  W,  and  Fe  [Figs.  1(a),  2(a),  and 
2(b)]  during  unloading  follow  almost  their  loading  curves. 
This  suggests  that  the  SE  of  Mo,  W  and  Fe  wires  is  charac¬ 
teristic  of  low-energy  dissipation,  which  is  quite  different 
from  the  SE  in  fee  nanowires. 11,12,14 

B.  Irreversible  deformation  in  V  nanowires 

Unlike  Mo,  W,  and  Fe  wires,  the  vanadium  (V)  nanowire 
failed  to  exhibit  any  SE  under  unixial  tension.  As  shown  in 
Fig.  2(c),  V  nanowires  first  deform  elastically  during  loading. 
Upon  yielding,  necking  occurs  and  this  quickly  leads  to  thin¬ 
ning  and  eventual  rupture  of  the  wire.  While  twinning/ 
detwinning  is  responsible  for  the  superelasticity  in  Mo,  W, 
and  Fe  wires,  (111)/{112}  dislocation  slip  is  primarily  re¬ 
sponsible  for  the  thinning  and  eventual  rupture  process  in  V 
nanowires  [shown  later  in  Fig.  5(d)].  The  slip  mechanism, 
involving  full  dislocations,  causes  the  tensile  deformation  to 
be  permanent  and  irreversible  upon  unloading.  Consequently, 
no  superelasticity  is  possible  for  V  nanowires. 

C.  Effect  of  strain  rate  on  the  superelasticity  in  Mo,  W,  and 
Fe  nano  wires 

Our  MD  simulations  are  performed  with  the  strain  rate  of 
107  s_1.  Although  it  is  already  very  slow  for  MD  simula¬ 
tions,  it  is  an  extremely  high  strain  rate  compared  to  reality. 
Previously,  work  related  to  fracture  have  reported  twinning 
as  a  mechanism  at  high  strain  rates  in  iron  but  with  slower 
rates  a  dislocation  mechanism  takes  over.26,27  To  verify 
whether  the  observed  superelasticity  mediated  by  the  twin¬ 
ning  mechanism  in  our  simulations  is  dependent  on  the  strain 
rate,  we  further  checked  the  deformation  behavior  of  Mo,  W, 
and  Fe  nanowires  under  quasistatic  conditions  by  using  an 
annealing  MD  method.  To  perform  this  calculation,  we  first 
applied  a  prescribed  strain  (0.014%)  to  the  relaxed  [100]- 
oriented  nanowires  through  the  fixed  atomic  layers  at  the 
ends  at  0  K  and  then  minimized  the  energy  of  the  nanowire 
using  the  conjugate-gradient  (CG)  method.  By  repeating  the 


FIG.  2.  (Color  online)  Stress-strain  curves  for  W,  Fe,  and  V 
nanowires  at  300  K.  Figures  2(a)  and  2(b)  show  the  superelasticity 
with  the  aid  of  twin  deformation  in  a  W  and  Fe  nanowire,  respec¬ 
tively.  Figure  2(c)  shows  the  failure  of  superelasticity  in  the  V  nano¬ 
wire  due  to  the  activation  of  full-dislocation  slip. 

above  process,  we  completed  the  tensile  deformation  of  the 
wires.  Unloading  follows  a  similar  path  by  applying  a  nega¬ 
tive  strain  increment. 

Figures  3(a)-3(c)  show  the  calculated  stress-strain  curves 
for  Mo,  W,  and_Fe  nanowires  which  have  dimensions  of 
50a X 5y2a X 5y2a,  where  a  is  lattice  constant  at  0  K.  It  is 
clear  that,  even  under  quasistatic  loading/unloading  at  0  K, 
W,  Mo,  and  Fe  nanowires  can  still  exhibit  almost  perfect  or 
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FIG.  3.  (Color  online)  Stress-strain _curves  of  Mo,  W,  and  Fe 
nanowires  with  dimensions  of  50aX5y2aX5y2a  using  annealing 
MD  methods,  (a)  and  (b)  show  that  Mo  and  W  are  able  to  exhibit 
good  superelasticity  with  more  than  40%  recoverable  strain,  and  (c) 
shows  Fe  nanowire  can  exhibit  partial  superelasticity. 

partial  SE.  Further,  atomic  configurations  show  that  the  SE 
of  W,  Mo,  and  Fe  nanowires  are  still  mediated  by  a  twinning/ 
detwinning  mechanism,  similar  to  the  cases  of  loading/ 
unloading  at  300  K  with  high  strain  rate  [as  shown  in  Fig. 
1(b)].  This  indicates  that  the  reversible  twinning  behavior  we 
found  in  W,  Mo,  and  Fe  nanowires  is  independent  of  the 
strain  rate.  In  addition,  we  noticed  that  the  recoverable  strain 


FIG.  4.  (Color  online)  Stress-strain  curves  for  a  (a)  5.3-nm- 
width  Mo  and  (b)  6.7-nm-width  W  nanowire  at  300  K.  Both  show 
excellent  superelasticity. 

in  W,  Mo,  and  Fe  nanowires  is  different  at  0  K,  indicating 
that  the  mobility  of  twins  in  these  three  metals  (Mo,  W,  and 
Fe)  at  nanoscale  is  diffident  and  we  will  discuss  it  later  (Secs. 
IV  B  and  IV  C). 

D.  Size  effect  on  the  superelasticity  in  Mo  and  W  nanowires 

The  present  calculations  already  show  SE  in  W,_Mo,  and 
Fe  nanowires  with  a  small  cross  sectional  size  (5V2 a,  about 
2.2  nm,  2.2  nm,  and  2.0  nm,  respectively).  It  is  interesting  to 
determine  whether  such  an  effect  can  be  found  in  larger 
sizes,  closer  to  realistic  cases.  For  this  purpose,  we  took  W 
and  Mo  nanowires  as  examples,  and  studied  the  effect  of  the 
cross-sectional  size  on  the_ superelasticity  at  300  K  by  keep¬ 
ing  the  aspect  ratio  as  5\2- 

Figure  4(a)  shows  that  the  stress-strain  curve  for  a  5.3- 
nm-width  Mo  nanowire  upon  loading  and  unloading.  We  ob¬ 
served  that  the  Mo  nanowire  still  exhibits  almost  perfect  su¬ 
perelasticity.  In  addition.  Fig.  4(b)  shows  that  a  W  nanowire 
can  also  exhibit  superelasticity  even  when  the  cross-sectional 
size  reached  6.7  nm.  As  the  current  ability  to  synthesize 
nanowires  with  tunable  diameters  has  reached  1.6  nm  (Ref. 
28)  or  even  down  to  approximately  0.8  nm,29  e.g.,  gold  nano¬ 
wires  formed  by  electron-beam  irradiation  in  an  ultrahigh 
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FIG.  5.  (Color  online)  The  atomic  images  of  the  incipient  plasticity  upon  loading  in  Mo,  W,  Fe,  and  V  nanowires.  Figure  5(a)-5(c)  show 
the  snapshots  of  twin  nucleation  in  Mo,  W,  and  Fe,  respectively  (only  two  adjacent  lattice  planes  are  shown  for  each  snapshot).  Using  Mo 
as  a  typical  example,  a  stacking  fault  forms  at  first  from  the  perfect  nano  wires.  The  embryo  subsequently  forms  in  a  layer-by-layer  manner 
with  respect  to  the  stacking  fault.  Yellow  lines  denote  the  stacking  layers  along  the  original  [100]  axis  of  the  nanowire.  Red  lines  refer  to  the 
newly  formed  microtwin.  Figure  5(d)  shows  the  full-dislocation  slip  in  the  V  nanowire,  leaving  a  step  at  the  surface  indicated  by  the  arrows. 


vacuum  electron  microscope  have  a  thickness  of  0.8-3  nm 
and  a  length  of  5-10  nm.29  The  present  size  (e.g.,  5.3  nm  or 
6.7  nm)  of  the  nanowire  already  reaches  dimensions  that  can 
be  perfectly  fabricated  in  experiments.28,29 

E.  Atomic  processes  for  twin  formation 
and  full  dislocation  slip 

From  the  above  discussions,  we  note  that  superelasticity 
in  Mo,  W,  and  Fe  nanowires  is  mediated  by  reversible  move¬ 
ment  of  twin  boundaries,  whereas  the  plastic  deformation  of 
V  nanowires  is  related  to  the  irreversible  slip  of  full  disloca¬ 
tions.  However,  we  would  like  to  understand  the  correspond¬ 
ing  atomic  processes  underlying  the  formation  and  growth  of 
twins  and  full  dislocation  slip.  We  therefore  consider  the  na¬ 
ture  of  the  atomic  configurations  at  the  initial  plastic  defor¬ 
mation  of  the  above  nano  wires. 

As  shown  in  Figs.  5(a)-5(c),  we  see  that  the  formation 
and  growth  of  twins  in  Mo,  W,  and  Fe  nanowires  is  quite 
similar:  the  nanowire  first  forms  a  one-layer  of  stacking 
faults  (SF)  from  the  perfect  crystal  through  a  |[111]  partial 
dislocation  in  the  {211}  plane,  followed  by  a  second  stacking 
fault  with  |[111]  partial  dislocation  in  the  adjacent  {211} 
plane.  These  layer-by-layer  stacking  processes  through  slid¬ 
ing  of  |[111]  partial  dislocations  in  adjacent  {211}  planes 
finally  cause  the  formation  and  growth  of  twins  in  Mo,  W, 
and  Fe  nanowires.  As  for  the  V  nanowire  [shown  in  Fig. 
5(d)],  we  see  that  the  first  SF  with  a  |[111]  partial  disloca¬ 
tion  is  also  easy  to  form  along  the  {21 1}  plane.  However,  the 
second  |[111]  partial  fails  to  be  activated  in  the  adjacent 


{211}  planes.  Rather,  another  f[lll]  partial  dislocation  is 
emitted  along  the  same  {211}  layer  and  results  in  the  forma¬ 
tion  of  a  full  dislocation  y[lll]  in  the  specific  {211}  plane. 

To  further  understand  the  growth  process  for 
(lll)/{112}-type  twins  in  bcc  nano  wires,  we  focus  on  the 
formation  of  SFs  during  tension  at  300  K.  The  results  show 
that  the  growth  of  <1 1 1)/{1 12}  twins  in  bcc  nanowires  in¬ 
volves  a  consecutive  movement  of  g[lll]  partial  disloca¬ 
tions  in  adjacent  {211}  planes.  As  illustrated  in  Fig.  6,  ini¬ 
tially  a  partial  dislocation  is  nucleated  at  one  side  of  the  free 
surface  of  the  Mo  nanowire  [Fig.  6(a)]  and  then  slides  along 
the  [111]  direction  in  the  {211}  planes  gradually  [Figs.  6(b) 
and  6(c)].  It  finally  moves  out  of  the  crystal  at  the  opposite 
side  leading  to  the  formation  of  a  one-layer  SF.  Subse¬ 
quently,  another  partial  dislocation  is  nucleated  in  a  similar 
way  [Fig.  6(d)]  to  generate  yet  another  partial  dislocation  in 
the  successive  {211}  plane.  Through  this  repeatable  process, 
the  growth  of  twinning  is  accomplished.  The  progress  of 
twin  motion  is  quite  similar  to  that  in  a  fee  nanowire,  which 
is  related  to  the  sliding  of  {[(112)  partial  dislocations  in  ad¬ 
jacent  {111}  planes.14 

IV.  DISCUSSION 

A.  Driving  force  for  superelasticity  of  bcc  nanowires 

It  is  well  known  that  twinning  in  bulk  materials  is  one  of 
the  deformation  modes  which  induces  permanent  plastic 
strain  upon  unloading.30,31  The  present  superelasticity  in  Mo, 
W,  and  Fe  nanowires  implies  that  detwinning  mediated  by 
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FIG.  7.  (Color  online)  The  size-dependent  surface  energy  differ¬ 
ence  per  atom  between  the  [100]-axis  and  [011]-axis  configurations 
in  Mo  nanowires  at  300  K.  With  decreasing  size,  the  energy  differ¬ 
ence  increases. 


(d) 40130  ps 


FIG.  6.  (Color  online)  Twinning  propagation  process  in  a  Mo 
nanowire  (only  two  adjacent  lattice  planes  are  shown  for  each  snap¬ 
shot).  A  |[111]  partial  dislocation  nucleates  on  one  side  of  the 
surface,  propagates  and  annihilates  in  the  (211)  planes.  The  con¬ 
secutive  partial  dislocation  at  free  surface  then  follows  a  similar 
process.  Red  arrow  indicates  the  slip  direction  of  the  partial  dislo¬ 
cation.  Green  line  denotes  the  twin  boundary.  Yellow  line  denotes 
the  symmetric  atomic  stacking  beside  the  twin  boundary.  Orange 
line  indicates  the  mismatch  of  atomic  layers  caused  by  the  move¬ 
ment  of  partial  dislocations. 


reverse  twin  boundary  motion  can  occur  spontaneously  dur¬ 
ing  unloading.  It  also  indicates  that  there  should  be  some 
intrinsic  force  to  drive  the  reversible  twin  boundary  motion. 

The  driving  force  for  spontaneous  reversible  twin  bound¬ 
ary  motion  in  Mo,  W,  and  Fe  nanowires  is  mainly  due  to  the 
reduction  of  surface  energy  of  the  nanowires,  similar  to  that 
in  fee  nanowires.11,12’32  Taking  Mo  as  an  example,  Fig.  7 


shows  the  averaged  surface  energy  per  atom  of  the  perfect 
[100]-axis  and  [011]-axis  Mo  nanowires  as  a  function  of  the 
cross-sectional  size  of  nanowire  at  300  K.  From  this  figure, 
we  clearly  see  that,  compared  to  the  original  [100]-axis  con¬ 
figuration,  the  reoriented  [011]-axis  configuration  has  a 
higher  surface  energy  per  atom,  providing  the  driving  force 
to  transform  back  to  the  [100]  configuration  upon  unloading. 
In  addition,  the  surface  energy  difference  between  the  two 
configurations  increases  nonlinearly  with  decreasing  cross 
section  of  the  nanowire.  This  indicates  that  superelasticity 
becomes  easier  to  induce  in  small  cross-section  nanowires 
and  the  shape  recovery  ability  becomes  weaker  or  even  van¬ 
ishes  for  larger  cross-section  nanowires. 


B.  Generalized  stacking  fault  energy  (GSFE)  of  twinning  and 
full  dislocation  slip  in  bcc  nanowires 

Our  simulation  results  have  shown  that  Mo,  W,  and  Fe 
nanowires  can  exhibit  SE  through  reversible  twinning, 
whereas  V  nanowires  do  not  show  such  behavior  due  to  the 
emission  of  full  dislocation  sliding.  Thus,  bcc  nanowires  that 
exhibit  SE  are  strongly  influenced  by  the  competition  be¬ 
tween  twinning  and  full  dislocation  slip.  Therefore,  it  is  natu¬ 
ral  to  ask  why  some  bcc  single-crystalline  nanowires  (e.g.. 
Mo,  W,  and  Fe)  can  undergo  twinning  and  not  the  full- 
dislocation  slip,  and  others  (e.g.,  V)  do  not  show  this  behav¬ 
ior. 

In  fee  nanowires  the  problem  can  be  understood  in  terms 
of  the  GSFE.33  36  However,  because  of  periodic  boundary 
conditions  in  directions  parallel  to  the  shear  planes,  the  in¬ 
fluence  of  surface  energy  change  during  the  shearing  process 
is  ignored.  As  a  result,  the  GSFE  derived  so  far  only  reflects 
the  characteristics  of  fee  bulk  materials  but  cannot  represent 
the  true  properties  of  nanoscale  systems.  Here  we  extend  this 
method  to  bcc  nanowires  by  including  the  effects  of  surface 
energy  changes  during  the  shearing  process  to  obtain  the 
GSFE  for  twinning  and  full-dislocation  slip.  We  first  estab¬ 
lish  a  model  for  twinning  and  full-dislocation  slip  in  bcc 
nanowires,  and  then  calculate  the  GSFE  of  twinning  and  full- 
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FIG.  8.  (Color  online)  A  sche¬ 
matic  illustration  of  full-dislocation 
slip  and  deformation  twinning  in  the 
(111)/{112}  system  of  bcc  nano¬ 
wires.  Atomic  configuration  of  a 
bcc  perfect  crystal  with  a  periodic 
stacking  sequence  as 

•  •  -ABCABCABC-  •  •  viewed  along 
[Oil]  directions.  Two  nearest  atoms 
in  the  horizontal  direction  have  a 
Burgers  vector  of  |[111].  (a)  Twin 
deformation.  For  each  step,  the  at¬ 
oms  in  the  red  part  shear  to  §[111], 
which  produces  the  same  displace¬ 
ment  as  in  the  adjacent  (211)  planes 
leading  to  the  stacking  sequence 
change,  (b)  Process  of  full  disloca¬ 
tion  slip  in  two  steps,  involving  the 
continued  emission  of  §[111]  and 
|[111]  partials  in  the  same  (211) 
plane.  After  sliding,  the  stacking  se¬ 
quence  has  no  change. 


dislocation  slip  in  Mo,  W,  Fe,  and  V  nanowires. 

To  begin  this  analysis,  we  first  confirm  the  deformation 
process  involving  the  two  competing  modes  (full-dislocation 
slip  and  twinning)  in  bcc  nanowires.  Based  on  the  simulation 
results  in  Fig.  5,  we  propose  a  model  to  explain  the 
(1 1 1)/{  1 12}-type  twinning  mechanism  of  bcc  nanowires,  a 
process  related  to  the  piling  of  §(111)  partial  dislocations. 
Here  we  just  use  a  specific  [111]/ (211)  system  to  illustrate 
these  processes.  The  stacking  sequence  of  (21 1)  planes  in  the 
bcc  system  can  be  regarded  as  •  •  -ABCABCABC-  ■  •  layers.  A 
full  dislocation  with  Burgers  vector  §[111]  can  dissociate 
into  two  partials 

ar  _  ar  _  ar 

-[111]  —  -[HI]  +  “[HI]. 

2  o  3 

The  nucleation  and  growth  of  twins  can  be  ascribed  to  the 
emission  of  §[111]  partial  dislocations  on  a  series  of  adja¬ 
cent  (211)  planes,  as  shown  in  Fig.  8(a)  schematically.  First, 
a  f[lll]  partial  slides  in  the  (211)  plane  to  produce  a  dis¬ 
placement  of  along  the  [111]  direction  to  form  a  stacking 
fault,  leading  to  the  alternating  stacking  sequence 

•  • -ABCBCABCA- • -.  Then  another  §[111]  partial  emits  in 
the  neighboring  (211)  planes  to  change  the  sequence  to 

•  •  -ABCBABCAB-  •  -.  The  subsequent  partials  similarly  fol¬ 
low  to  further  allow  for  twin  nucleation  and  growth.  In  the 


present  work,  we  calculated  the  GSFE  curve  for  twinning  up 
to  three  layers. 

Next,  we  determine  the  progress  of  full  dislocation  slid¬ 
ing.  It  is  known  that  there  are  three  main  dislocation 
slip  systems  in  bcc  metals,  i.e.,  (111)/{112}  [specifically,  a 
§[111]  full  dislocation  slip  along  the  (211)  plane], 
(111)/{110}  and  (111)/{123}. 37,38  Since  all  the  sliding  pro¬ 
cesses  follow  similarly,  here  we  only  illustrate  with  the 
(111)/{112}  slip.  As  shown  in  Fig.  8(b),  the  emission  of  a 
§[111]  full  dislocation  can  be  regarded  as  leading  §[1 1 1]  and 
trailing  §[111]  partial  slip  in  the  same  (211)  plane.  After 
emitting  the  full  dislocation,  the  atomic  stacking  sequence 
does  not  change;  it  still  remains  the  perfect 
•  •  -ABCABCABC-  ■  •  sequence. 

On  the  basis  of  the  above  models  of  twinning  and  full- 
dislocation  slip,  we  can  calculate  the  GSFE  for  bcc  nano¬ 
wires.  We  first  create  the  bcc  nanowires  with  50a  X  5y/2a 
X5y2 a  (a  is  the  lattice  constant  at  0  K),  the  same  dimen¬ 
sions  as  the  initial  model  of  the  present  MD  simulations. 
Shear  is  then  applied  to  the  nanowires  to  reproduce  the 
§[111]  full  dislocation  slip  along  {112},  {110},  and  {123} 
planes  and  the  (111)/{112}  twin  according  to  the  dynamical 
process  shown  in  Figs.  8(a)  and  8(b).  During  this  process, 
atoms  in  the  two  directions  that  are  perpendicular  to  the 
shear  direction  are  relaxed  by  means  of  CG  algorithm.  As  a 
result,  we  obtain  the  changes  in  potential  energy  with  the 
shear  displacement  r.  Taking  the  (111)/{112}  full  dislocation 
slip  and  twin  of  Mo  nanowire  as  an  example  (shown  in  Fig. 
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FIG.  9.  (Color  online)  Potential  energy  change  with  respect  to 
the  number  of  111]  partials  between  twining  (black  solid  line) 
and  full-dislocation  slip  (red  solid  line)  for  (J_l  l)/{ l_l 2}  systems  of 
Mo  nanowires  with  dimensions  of  50a  X  5  V2a  X  5  v2 a,  where  a  is 
the  lattice  constant  at  0  K.  The  dashed  lines  are  the  reference  for 
eliminating  the  tilting  of  the  curve.  The  blue  dashed  line  is  the 
reference  line  to  emit  the  first  partial  as  well  as  for  full-dislocation 
slip.  The  green  dashed  line  is  the  reference  line  for  the  emission  of 
the  second  and  the  third  partials  for  twinning. 

9),  we  note  that  the  calculated  potential  energy-shear  dis¬ 
placement  curves  are  no  longer  periodic  as  that  of  bulk 
materials, 33,39-41  instead  they  are  tilted  toward  a  high-energy 
state  due  to  the  newly  created  ledge  in  the  side  surface  of  the 
wire.  To  obtain  the  energy  change  due  to  the  formation  of  a 
stacking  fault,  one  needs  to  eliminate  the  energy  induced  by 
the  newly  created  ledge  in  the  surface.  Here,  we  used  a  simi¬ 
lar  method  as  in  Ref.  42,  to  eliminate  such  effect  by  subtract¬ 
ing  the  inclined  lines  (blue  dashed  line  in  Fig.  9)  from  the 
original  tilted  curves.  For  the  (111)/{112}  twin  (black  solid 
line  in  Fig.  9),  the  energy  change  due  to  the  first  ^[1 11] 
partial  should  be  the  same  as  the  (11 1)/{1 12}  full  dislocation 
(red  solid  line  in  Fig.  9)  slip  because  they  follow  the  same 
dynamical  process.  For  the  second  and  third  |[111]  partial 
dislocations,  an  inclined  line  that  connects  the  starting  and 
finishing  points  (green  dashed  line  in  Fig.  9)  is  also  used  as  a 
reference  to  remove  the  effect  of  newly  created  ledges.  The 
obtained  curves  after  the  above  treatment  were  then  normal¬ 
ized  by  dividing  the  energy  by  their  original  sliding  areas  in 
the  nanowire.  Finally  we  plot  the  GSFE  curves  with  an  inte¬ 
ger  number  of  g[lll]  partial  distance. 

Figures  10(a)-10(d)  show  the  calculated  GSFE  curves 
with  the  integer  number  of  f[lll]  partial  distance  for  Mo, 
W,  Fe,  and  V  nanowires.  The  curves  in  different  crystals 
show  similar  morphology.  Using  the  Mo  nanowire  as  an  ex¬ 
ample,  the  single  solid  line  and  three  dashed  lines  in  Fig. 
10(a)  represent  twinning  deformation  and  three  full- 
dislocation  slip  systems,  respectively.  We  see  that  the  GSFE 
curve  of  the  (111)/{112}  twin  can  be  divided  into  two  stages: 
First,  there  is  a  high  energy  barrier  ( ysf)  for  emitting  the  first 
|[111]  partial  to  form  a  SF  from  the  initial  perfect  crystal. 
The  shape  of  the  GSFE  curve  for  the  emission  of  subsequent 
|[111]  partials  is  similar  to  a  periodic  wave  with  peaks  and 
valleys.  We  define  the  difference  between  a  peak  and  a  valley 
of  the  wave  as  Ay""1.  Because  Aywin  is  the  energy  barrier  of 
emission  of  subsequent  |[111]  partials  on  the  adjacent  (211) 
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planes  in  the  presence  of  a  pre-exiting  SF,  it  is  related  to  the 
energy  barrier  for  the  subsequent  formation  of  twin  layers. 
From  Fig.  10(a),  we  also  see  that  the  GSFE  curves  among 
different  dislocation-slip  systems  are  similar,  there  is  a  peak 
in  the  GSFE  curve  indicating  the  maximum  energy  barrier 
for  the  emission  of  a  full  dislocation  from  the  perfect  crystal, 
and  we  define  the  peak  value  as  the  unstable  stacking  fault 
energy  -ys.  Besides,  the  difference  in  y"s  and  yf  indicates  the 
energy  barrier  to  form  a  full  dislocation  in  the  presence  of  a 
pre-exiting  SF,  and  we  define  it  as  Aysllp.  In  addition,  we 
note  that  the  GSFE  curve  for  full  dislocation  slip  in  the 
(111)/{112}  system  is  the  same  as  that  for  (111)/{112}  twin¬ 
ning  before  the  emission  of  the  first  ^[111]  partial  disloca¬ 
tion. 

C.  Competition  between  twinning  and  full-dislocation  slip  in 
bcc  nanowires 

Based  on  the  above  GSFE  curves,  we  are  now  ready  to 
investigate  the  competition  between  twinning  and  full- 
dislocation  slip  in  bcc  nanowires  and  understand  the  differ¬ 
ence  in  deformation  in  Mo,  W,  Fe,  and  V  nanowires.  Our 
analyses  show  that  a  SF  generated  by  the  |[111]/(211)  par¬ 
tial  dislocation  is  always  preferred  over  the  (111)/{110}  and 
(111)/{123}  full-dislocation  slip  systems  in  bcc  nanowires. 
However,  whether  (111)/{112}  twinning  or  ( 1 1 1 ) / { 1 1 2}  full- 
dislocation  occurs  in  bcc  nanowires  depends  on  the  compe¬ 
tition  between  subsequent  emission  of  7 [111]  partials  in  the 
neighboring  (211)  planes  and  the  emission  of  f[lll]  in  the 
same  (211)  plane.  If  the  former  is  selected,  then  (111)/{11 2} 
twin  is  preferred;  otherwise  (111)/{112}  full  dislocation  slip 
will  be  stimulated  in  the  bcc  nanowires. 

We  first  discuss  why  Mo,  W,  and  Fe  bcc  single-crystalline 
nanowires  can  undergo  SE  mediated  by  twinning  deforma¬ 
tion  whereas  V  nanowires  can  only  endure  permanent  plastic 
deformation  through  full-dislocation  slip.  According  to  our 
previous  analysis,  the  necessary  condition  to  form  a  twin 
from  the  perfect  nanowire  is  the  formation  of  stacking  fault 
emitted  by  one  |[111]/(211)  partial  dislocation.  This  leads 
to  a  tendency  to  form  a  twin  by  further  emission  of  |[111] 
partial  dislocations  in  adjacent  (211)  planes.  On  the  other 
hand,  if  one  of  the  three  full-dislocation  slip  systems  is  acti¬ 
vated  from  the  perfect  crystal,  the  bcc  perfect  crystal  should 
overcome  the  unstable  stacking  fault  (y*s).  Obviously,  the 
slip  system  with  min{ys},  which  is  a  minimum  of  y“s  among 
the  ysl10,  ysl12,  ysl23  in  <111)/{110},  (11 1)/{1 12},  and 
(11 1) /{ 123}  slip  systems,  is  more  likely  to  be  activated. 
Since  yf  and  min{y*s}  characterize  the  difficulty  of  forming 
twinning  partial  dislocation  and  full-dislocation  slip,  we  then 
use  the  ysi/niin{yus}  ratio  to  describe  the  competition  be¬ 
tween  these  two  deformation  modes.  Clearly,  if 
y^/minly5}  >  1 ,  full-dislocation  slip  is  preferred;  if 
ysf/min{yus}<  1,  a  stacking  fault  induced  by  a  |[111]/(211) 
partial  dislocation  is  easier  at  the  initial  stage  of  deformation. 

Table  I  lists  the  calculated  y“s  for  different  slip  systems, 
ysf,  A/Wln,  Aysllp,  and  ysf/ min{  y“s}  for  W,  Mo,  Fe,  and  V 
nanowires.  To  make  a  comparison,  we  also  calculated  the 
above  parameters  in  bulk,  which  are  based  on  the  traditional 
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FIG.  10.  (Color  online)  GSFE  of  BCC  Mo,  W,  Fe,  and  V  crystals,  respectively,  shown  in  (a)-(d).  The  solid  and  dashed  lines  represent 
twinning  and  full-dislocation  slip,  respectively,  yf  refers  to  the  barrier  for  producing  stacking  fault  energy  by  emission  of  the  first  f[lll] 
partial,  y“s  represents  the  unstable  stacking  fault  energy  to  induce  the  full-dislocation  slip,  which  is  specifically  defined  as  y^110,  y^112,  and 
yisi23  jn  ^  jjjj-gg  (1H)/{110},  (111)/{112},  and  (111)/{123}  systems,  respectively.  Ay1*111  and  Aysllp  refer  to  the  maximum  energy  barrier 
by  emission  of  the  subsequent  partials  for  twin  growth  and  full  dislocation  slip  in  the  presence  of  pre-existing  stacking  fault  of  the 
(111)/{112}  system. 


methods,41  i.e.,  using  the  periodic  boundary  conditions  along 
the  directions  parallel  to  the  shear  plane.  As  shown  in  Table 
I,  we  observe  that  the  difference  between  ys  in  (111)/{1 10}, 
(111)/{112},  and  (111)/ {123}  slip  systems  and  yf  for  nano¬ 
wires,  and  in  bulk  is  small  in  W,  Mo,  and  Fe,  but  a  little 
larger  in  V.  This  results  in  the  changes  in  ysl/ min{ylls}  and 
Aytwln/Ayslip  in  W,  Mo,  Fe,  and  V  nanowires.  This  is  reason¬ 
able  as  the  surface  energies  play  a  crucial  role  at  the 
nanoscale. 5-9,42 

As  the  present  paper  focuses  on  the  deformation  behavior 
at  nanoscale,  we  use  the  data  obtained  from  the  nanowires  to 
analyze  the  competition  between  full-dislocation  slip  and 
twinning.  As  shown  in  Table  I,  the  values  of  yif/min{yus}  in 
Mo,  W,  Fe,  and  V  nanowires  are  all  smaller  than  1  (0.73, 
0.79,  0.98,  and  0.93,  respectively).  This  indicates  that  in  bcc 
nanowires  a  §[111]  partial  dislocation  is  preferred  over  the 
(111)/{110}  and  <111)/ {123}  full  dislocation  slip  at  the  ini¬ 
tial  stage  of  deformation. 

After  the  SF  induced  by  one  g[lll]/(211)  partial  dislo¬ 
cation  is  formed,  it  is  not  clear  whether  <11 1)/{1 12}  twin  or 
full  dislocation  will  occur  from  the  perfect  crystal.  It  depends 
on  the  competition  between  the  energy  barrier  to  emit  a 
|[111]  trailing  dislocation  (Aysllp)  to  form  a  full  dislocation 


and  the  energy  to  emit  a  subsequent  |[111]  partial  disloca¬ 
tion  (Aytwln)  to  form  a  twin.  Clearly,  the  ratio  of 
Aytwln/Aysllp  reflects  a  competition  between  twinning  and 
full-dislocation  slip,  if  Aywin/ Aysllp>  1,  it  means  the  emis¬ 
sion  of  a  f[lll]  trailing  dislocation  on  the  pre-existing  SF  is 
easier,  and  thus  <111)/{1 12}  full  dislocation  is  preferred; 
whereas  if  Aytwm/Aysllp  <  1 ,  the  emission  of  a  subsequent 
^[111]  partial  dislocation  on  the  pre-existing  SF  is  easier, 
and  then  the  <111)/{1 12}  twin  is  preferred. 

As  shown  in  Table  I,  the  values  of  Aywin/Aysllp  for  Mo, 
W,  and  Fe  nanowires  are  all  smaller  than  1  (0.06,  0.23,  and 
0.42,  respectively).  This  implies  that  the  <11 1)/{1 12}  twin  is 
preferred  in  Mo,  W,  and  Fe  nanowires.  However,  the  value  of 
Aywm/Ayllp  in  V  nanowires  is  larger  than  1  (1.15),  which 
indicates  that  the  |[111]/(211)  full  dislocation  is  preferred 
in  V  nanowires.  This  analysis  is  also  consistent  with  the 
atomic  configurations  in  Figs.  5(a)-5(d). 

Based  on  the  above  understanding,  we  can  further  predict 
the  preference  for  twinning  in  bcc  crystals.  As  twinning  is 
related  to  the  formation  of  the  SF  by  a  |[111]/(211)  partial 
dislocation  and  then  subsequent  formation  of  |[111]  partials 
on  adjacent  (211)  planes,  it  is  natural  to  conclude  that  due  to 
a  smaller  value  of  yfyminly18}  and  Aywin/Aysllp,  it  is  easier 
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TABLE  I.  Calculated  ys  for  different  slip  systems,  yf,  Ay™,  Aylip,  and  yf/  min{yl‘'}  in  bulk  and  in 
nanowires  (with  size  50a  X  5\l2a  X  5\2a,  where  a  is  the  lattice  constant  at  0  K)  of  W,  Mo,  Fe,  and  V.  The 
energy  unit  is  mj/m2. 


w 

Mo 

Fe 

V 

Nano 

Bulk 

Nano 

Bulk 

Nano 

Bulk 

Nano 

Bulk 

yasllO 

1057 

1021 

952 

902 

625 

660 

965 

1075 

yiS112 

1144 

1185 

1038 

1063 

738 

766 

1069 

1238 

yaS123 

1083 

1111 

990 

1043 

679 

754 

973 

1243 

/ 

773 

830 

754 

796 

615 

649 

900 

960 

Ay™ 

21 

61 

66 

96 

52 

63 

195 

198 

Ay  ^ 

371 

355 

284 

266 

123 

117 

169 

278 

yf/min(/s) 

0.73 

0.81 

0.79 

0.88 

0.98 

0.98 

0.93 

0.89 

Ay™/Aylip 

0.06 

0.17 

0.23 

0.36 

0.42 

0.54 

1.15 

0.71 

for  the  system  to  form  a  (111)/{112}  twin.  As  shown  in 
Table  I,  both  the  ratio  of  y^/minjy15}  and  Aywl7Aysllp  in¬ 
dicate  the  sequence  W  <Mo<Fe<  1.  Thus,  compared  to  full 
dislocation  slip  systems,  the  preference  for  twinning  in  the 
present  work  is  given  by  W  >Mo>Fe. 

D.  Comparison  of  superelasticity  between  fee  and 
bcc  nanowires 

We  have  demonstrated  above  that,  driven  by  surface  en¬ 
ergy  minimization,  bcc  nanowires  can  show  superelasticity 
through  twinning  deformation.  Combined  with  the  work  for 
fee  systems, 10-12  our  present  study  shows  that  SE  should  be  a 
general  phenomenon  in  cubic  metal  nanowires  under  the 
condition  of  high  surface  energy  and  twinning  deformation. 
The  minimization  of  surface  energy  provides  an  intrinsic 
driving  force  for  the  transformation  between  two  configura¬ 
tions,  and  twinning  deformation  provides  the  dynamic  path¬ 
way  for  the  completion  of  crystal  reorientation.  However,  it 
is  interesting  to  examine  the  differences  between  bcc  and  fee 
nanowires  in  their  superelastic  behavior.  Here  in  points  (1)- 
(3)  below,  we  show  that  although  bcc  nanowires  have  similar 
shape  recovery  behavior,  they  behave  differently  from  their 
fee  counterparts  regarding  twin  formation,  twin  mobility,  en¬ 
ergy  dissipation,  and  applications  over  a  wider  temperature 
range. 

(1)  The  GSFE  for  twinning  in  bcc  metals  is  obviously 
different  from  that  of  fee  metals,  indicating  a  dramatic  dif¬ 
ference  in  twinning  formation  and  mobility  between  bcc  and 
fee  metals.  First,  the  energy  barrier  for  emitting  the  first  par¬ 
tial  to  form  a  SF  in  a  bcc  metal  is  quite  large  compared  to 
that  in  a  fee  metal.  For  example,  the  barrier  is  773  mJ/m2 
and  180  mJ/m2  in  W  and  Cu,14  respectively.  Since  the  for¬ 
mation  of  SF  is  of  great  importance  for  nucleation  of  twins, 
bcc  metals  need  a  large  driving  force  for  twinning  nucleation 
compared  to  fee  metals.  Second,  the  energy  barrier  for  the 
emission  of  subsequent  partial  dislocations  for  twinning  in 
bcc  metals  is  smaller  than  that  for  a  full  dislocation.  In  con¬ 
trast,  the  energy  barrier  for  the  emission  of  subsequent  partial 
dislocations  for  twining  in  fee  metals  is  larger  than  that  for  a 
full  dislocation.  '4  Therefore,  in  bcc  metals  there  is  a  bigger 


opportunity  for  the  emission  of  subsequent  partial  disloca¬ 
tions,  making  the  formation  of  twins  much  easier  than  that  in 
fee  metals.  Then  the  shape  of  GSFE  curves  after  SF  forma¬ 
tion  in  bcc  metals  is  also  significantly  different  from  that  of 
fee  metals.  The  energy  landscape  for  fee  twinning  is  always 
very  rugged;  in  contrast  to  the  flat  energy  landscape  in  bcc 
metals.  Thus,  bcc  metals  have  a  very  low  twin  boundary 
migration  energy  and  twin  boundary  motion  in  bcc  metals  is 
easier  than  in  fee  metals  [e.g.,  21  mJ/m2  in  W  vs 
157  mJ/m2  in  Cu  (Ref.  14)]. 

(2)  The  difference  in  twinning  mobility  results  in  obvious 
differences  in  stress-strain  curves  between  bcc  and  fee  met¬ 
als.  In  contrast  to  the  typical  strain  hardening  behavior  in  fee 
nanowires,  there  is  no  obvious  strain  hardening  in  the  present 
study  of  bcc  systems.  Most  interestingly,  the  stress-strain 
curves  of  the  bcc  nanowires  show  quite  a  small  hysteresis 
and  lower  energy  dissipation  during  the  loading/unloading 
process  that  is  significantly  different  from  that  in  fee  crystals. 
This  difference  is  easily  understood  in  terms  of  the  difference 
of  twin  boundary  motion  barrier  between  bcc  and  fee  metals. 
Due  to  the  very  low  twin  boundary  barrier,  the  twins  in  bcc 
nanowires  move  easily  and  revert  back  than  their  fee  coun¬ 
terparts,  leading  to  the  small  strain  hardening  rate  and  small 
hysteresis. 

(3)  Some  refractory  bcc  nanowires,  such  as  W  and  Mo, 
can  show  SE  in  quite  a  wide  temperature  range  than  that  in 
fee  nanowires,  fee  nanowires,  e.g.,  Cu,  can  only  show  supe- 
realsticity  at  temperatures  lower  than  900  K  due  to  the  low 
melting  point.11  Above  this  temperature,  surface  diffusion 
becomes  intense,  and  destroys  the  nanostructures.  However, 
as  bcc  W  and  Mo  crystals  have  extremely  high  melting 
points  (3695  K  and  2896  K,  respectively),  surface  diffusion 
is  not  that  important  at  relatively  high  temperature  for  these 
refractory  metals.  We  find  that  W  nanowires  can  show  SE 
even  at  temperatures  higher  than  1500  K  (shown  in  Fig.  11), 
which  exceeds  the  value  for  almost  all  the  reported  high- 
temperature  shape  memory  alloys  in  bulk. 16,17  Such  behavior 
may  help  the  development  of  high-temperature  shape- 
memory  materials  at  the  nanoscale. 

Finally,  it  should  be  noted  that,  our  further  work43,44  will 
show  that,  some  hexagonal  closed  packed  structures,  e.g.,  Zr 
metal  nanowires,  can  also  show  superelasticity  driven  by  the 
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FIG.  1 1 .  (Color_online)_Superelasticity  in  a  W  nanowire  of  di¬ 
mensions  50fl X 5<2a X 5v2a  at  high  temperature  (1500  K),  where 
a  is  the  lattice  constant  at  0  K. 

reduction  of  surface  energy  but  mediated  by  an  inverse  mar¬ 
tensitic  transformation.  This  is  a  pathway  for  surface  energy 
minimization  that  completely  differs  from  that  in  fee  and  bcc 
metal  nanowires.  Combined  with  the  present  work  on  bcc 
nanowires  and  previous  studies  on  fee  nanowires,  it  clearly 
shows  that  superelasticity  should  be  even  more  general  in 
metal  nanowires  and  the  driving  force  comes  from  the  same 
origin — the  reduction  in  surface  energy  of  the  nanowire. 


and  twinning  deformation.  The  minimization  of  surface  en¬ 
ergy  provides  an  intrinsic  driving  force  for  the  transforma¬ 
tion  between  two  configurations  and  twinning  deformation 
provides  the  dynamical  path  for  the  complete  crystal  reori¬ 
entation.  (2)  In  bcc  nanowires,  the  (111)/{112}  stacking  fault 
generated  by  7(111)  partial  dislocations  is  always  preferred 
over  the  (111)/{110}  and  <111)/ {123}  full  dislocation  slip. 
However,  whether  a  <11 1)/{1 12}  twin  or  a  <11 1)/{1 12}  full 
dislocation  slip  is  preferred  in  bcc  nanowires  depends  on  the 
competition  between  the  emission  of  a  subsequent  §(111) 
partial  dislocation  on  adjacent  {112}  planes  and  the  emission 
of  a  |(1 1 1)  partial  in  the  same  plane.  If  the  energy  barrier  for 
emission  of  consecutive  |(111)  partials  adjacent  to  the  pre¬ 
existing  stacking  fault  is  lower  than  that  for  a  f(lll)  partial 
on  the  pre-existing  SF,  a  (11 1)/{1 12}  twin  is  preferred.  (3) 
The  generalized  stacking  fault  energy  for  twinning  in  bcc 
metals  is  quite  different  than  that  of  fee  metals.  Thus,  the  bcc 
metal  has  a  higher  energy  barrier  for  the  nucleation  of  a  twin, 
but  a  lower  energy  barrier  for  twin  formation  and  migration, 
which  results  in  some  unique  characteristics  of  the  superelas¬ 
ticity  in  bcc  nanowires,  such  as  low-energy  dissipation  and 
low  strain  hardening.  (4)  Due  to  their  extremely  high  melting 
points,  refractory  Mo  and  W  nanowires  can  show  superelas¬ 
ticity  at  very  high  temperatures,  exceeding  the  temperature 
for  fee  shape-memory  nanowires  and  almost  all  the  reported 
high-temperature  shape  memory  alloys  in  bulk. 


V.  CONCLUSION 

We  have  used  molecular-dynamics  simulations  to  study 
the  mechanical  behavior  of  bcc  single-crystalline  nanowires. 
Our  main  conclusions  can  be  summarized  as  follows:  (1) 
driven  by  surface  energy  minimization,  bcc  nanowires  (Mo, 
W,  and  Fe)  can  show  superelasticity  through  twin  deforma¬ 
tion.  The  superelasticity  in  bcc  nanowires,  in  association 
with  the  observed  superealsticity  in  fee  nanowires, 10-12  sug¬ 
gests  that  superelasticity  is  a  general  phenomenon  in  cubic 
metal  nanowires  under  the  conditions  of  high  surface  energy 
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